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Abstract 

We derive the equations of motion of spinning compact binaries including the spin-orbit (SO) 
coupling terms one post-Newtonian (PN) order beyond the leading-order effect. For black holes 
maximally spinning this corresponds to 2.5PN order. Our result for the equations of motion 
essentially confirms the previous result by Tagoshi, Ohashi and Owen. We also compute the spin- 
orbit effects up to 2.5PN order in the conserved (Noetherian) integrals of motion, namely the 
energy, the total angular momentum, the linear momentum and the center-of-mass integral. We 
obtain the spin precession equations at 1PN order beyond the leading term, as well. Those results 
will be used in a future paper to derive the time evolution of the binary orbital phase, providing 
^ \ more accurate templates for LIGO-Virgo-LISA type interferometric detectors. 
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I. INTRODUCTION 



The laser interferometer gravitational-wave (GW) detectors LIGO, Virgo, GEO 600 and 
TAMA300 are currently searching for GWs emitted by inspiralling compact binaries com- 
posed of neutron stars and/or black holes. Analyzing the data usingmatched filtering tech- 
nique requires a high-precision modelling of the inspiral waveform [l|, 0, 0, 0, H, @, 0, B 0] . 
The post-Newtonian (PN) approximation to general relativity has been applied to build 
accurate theoretical templates up to the 3.5PN precision level 1 for non-spinning compact 



bodies [10|, QjJ, [12|. Post-Newtonian templates are currently used in analysing the data 
with ground-based detectors and in the future they will be used to detect GWs emitted by 
supermassive black-hole binaries with the space-based detector LISA. 

Astrophysical observations suggest that black holes can have non- negligible spins, e.g., 
due to spin up driven by accretion from a companion during some earlier phase of the 
binary evolution. For a few black holes surrounded by matter, observations indicate a 



significant intrinsic angular momentum (see, e.g., Refs. [13j, [14], |l5j for stellar black holes 
and Refs. [TJ for supermassive black holes). The spin may even be close to its maximal 
value 18fl. Very little is known however about the black-hole spin magnitudes in binary 
systems [191 ] . 

To successfully detect GWs emitted by spinning, precessing binaries and to estimate the 
binary parameters, spin effects should be included in the templates. For maximally spinning 
compact bodies the spin-orbit coupling (linear in the spins) appears dominantly at the 1.5PN 
order, while the spin-spin one (which is quadratic) appears at 2PN order. The spin effect on 
the free motion of a test particle was first obtained in the form of a coupling to curvature 

Seminal works by Barker and O'Connell (23l . 24j yielded 
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by Papapetrou et al. [2 
both the leading order spin-orbit and spin-spin contributions in the PN equations of motion. 

leading spin-orbit and spin-spin 
261 and predictions for spin-spin 



Based on the works 
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More recently, using an effective field theory approach [25 
couplings in the two-body Hamiltonian were re-derived 
couplings at 3PN order in the spin potential were obtained 

Kidder, Will and Wiseman [23, H^] (see also Refs. [3(J |31() computed the corresponding 
coupling terms in the radiation field, enabling thereby the derivation of the orbital phase 
evolution, the latter being the crucial quantity that determines the templates. Currently, 
only the leading order spin effects, i.e., the spin-orbit and spin-spin couplings have been 
implemented in the templates for spinning, precessing black-hole binaries [8|, [32|, 



371 ] started the computation of the 1PN 
linear in the spins, ap- 



More recently, Tagoshi, Ohashi and Owen [3| 
corrections to the leading spin-orbit coupling. Those corrections, 
pear at 2.5PN order. However, their work has never been completed: the very important 
conserved integrals associated to the equations of motion at 2.5PN order and the mass 
quadrupole moment at the 2.5PN order were not computed. 

The aim of the present paper together with its companion 38j is to complete the work 
of Refs. I 



371 ] and get the orbital phase evolution at 2.5PN order. In this paper we derive 
confirming the main result of Ref. 



the equations of motion, confirming the main result of Ref. 37j (but correcting several 
important misprints) and compute the entire set of conserved Noetherian integrals of the 
motion associated with the Poincare invariance, notably the energy and the total angular 
momentum. In Ref. [38[ (henceforth paper II) we evaluate the multipole moments and the 



1 As usual nPN refers to terms of order (v/c) 2n where v is the internal velocity and c the speed of light. In 
this paper we explicitely display all powers of c and of Newton's constant G. 
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radiation field so as to deduce the orbital phase evolution. 

The spin of a rotating body is of the order S true ~ mav sp i n , where m and a denote the 
mass and typical size of the body respectively, and where w sp i n represents the velocity of the 
body's surface. Here, by "true", we mean that the spin we are referring to is not rescaled 
[as in Eq. (11.11) below]. In this paper we shall consider bodies which are both compact, 
a ~ and maximally rotating, t> spin ~ c. For such objects the magnitude of the spin is 

roug hly S tTne The previous estimate shows that the spin goes as one power of 1/c, 

i.e., from the PN point of view, it is formally of order 0.5PN. Again, such a counting is 
appropriate for maximally rotating compact objects. It is then also customary to introduce 
a dimensionless spin parameter, generally denoted by Xi defined by S true = ^-^X- m our 
computation the use of such parameter x is n °t very convenient because it forces us to 
introduce some unwanted powers of the mass in front of the spins. On the other hand, it is 
useful to keep track of the correct PN order by counting all the powers of 1/c. Accordingly 
we shall "artificially" make explicit the factor 1/c in front of the spin by posing S trnc = S/c 
where S will be considered to be of "Newtonian" order. Hence, we shall denote the spin 
variable by 

S = c S tmc = Gm 2 X - (1-1) 

Such a notation displays explicitly all powers of 1/c for maximally rotating compact objects. 
Notably, the spin-orbit (SO) effect always carries a factor 1/c 3 in front, so that it is regarded 
as being of order 1.5PN, while the spin-spin (SS) effect appears at order 2PN in our termi- 
nology, and the 1PN correction to the spin-orbit is 2.5PN order. This PN counting for spin 
effects corresponds to the standard practice when defining the templates of LIGO/Virgo and 
LISA detectors (see Refs. 0,11]). 

For slowly rotating compact objects (f sp m ^ c ) the spin is formally of higher order, 
namely S true ~ Gm J s P in ^ 1/c 2 , hence the spin-orbit and spin-spin couplings are pushed at 
the 2PN and 3PN levels respectively. The 1PN correction to the spin-orbit manifests itself 
at the same level as the spin-spin coupling, namely 3PN. Of course all the computations 
in this paper and paper II [38j] are still valid in the case of slow rotation, but in this case 
the spin terms are expected to be numerically smaller, and comparable to higher-order PN 
contributions. 

This paper is organized as follows. In Sec. [TT] we describe the stress-energy tensor of 
spinning point particles and review some relevant features of the spin formalism before 
defining our spin variables. In Sec. II III we recall some general expressions of the PN metric 
and equations of motion, which are valid for arbitrary extended matter configurations. The 
PN metric is parametrized by certain elementary potentials computed in Sec. IIVI Our 
final results for the spin-orbit terms in the equations of motion at the 2.5PN order are 
presented in Sec. |V] in a general frame. They are also specialized to the center-of-mass 
frame and reduced to circular orbits. The precessional equations for the spins including the 
1PN relative correction are derived in Sec. I VI I Finally, in Sec. IVII[ we obtain the spin- 
orbit contributions to the conserved integrals associated with our 2.5PN dynamics. The two 
Appendices are devoted to some tests of our results. 



II. STRESS-ENERGY TENSOR FOR SPINNING POINT-PARTICLES 



Our calculations are based on the standard model of point-particles with spins |20|, [21 
21 0, 0, [H, 0, El, H 0, H H S]- In the Dixon formulation 0, the stress-energy 
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tensor, 

M S 

is the sum of the "monopolar" (M) piece, which is a linear combination of monopole sources, 
i.e. made of Dirac delta-functions, plus the "dipolar" or spin (S) piece, made of gradients 
of Dirac delta-functions. The four- dimensional formulation of the monopolar part reads as 

^^Ef^ ^P 1 , (2-2) 

M A J-oo V~9A 

where 5^' is the four- dimensional Dirac function. The world-line of particle A (A=l,2), 
denoted y A , is parametrized by the particle's proper time ta- The four-velocity is given by 
cu A = dy A /dr A , and normalized to g^ u u A u A = — 1, where g^ u = g^iljA) denotes the metric 
at the particle's location (the determinant of the metric at point A being denoted by g A )- 
The four-vector p A is the particle's linear momentum satisfying Eqs. (I2.4p — (j2.5j) below. The 
dipolar or spin part of the stress-energy tensor, which vanishes in the absence of spins, is 2 



b A 



-oo 



dr A S^u ^ ){X -l A) 
y—QA 



(2.3) 



where V p is the covariant derivative associated with the metric g^ at the field point x, and 
the anti-symmetric tensor S 1 ^ represents the spin angular momentum for particle A. 
The momentum-like quantity p A is a time-like solution of the equation 

= cu&pSZ = c 2 {f A u A - P X) . (2.4) 

The equation of motion of the particle with spin, equivalent to the covariant conservation law 
of the total stress-energy tensor, namely V V T^ U = 0, is given by the Papapetrou equation [2p| . 
22, Q 



Df A 1 



A Q^P„.v 



dTA 2 S>^ AMp . (2.5) 

The Riemann tensor is evaluated at the particle's position A, R AuXp = R^xpiVA)- The 
equation of motion (12.51) can also be derived directly from the action principle of Bailey and 
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It is well-known that a choice must be made for a supplementary spin condition (SSC) in 
order to fix unphysical degrees of freedom associated with some arbitrariness in the definition 
of S^ v . This arbitrariness can be interpreted, in the case of extended bodies, as a freedom in 
the choice for the location of the center-of-mass worldline of the body, with respect to which 
the angular momentum is defined (see e.g. [2t| for discussion). In this paper we adopt the 
covariant supplementary spin condition 

STPu=0, (2-6) 



2 Recall that with our convention the spin variable has the dimension of a true spin times c; the stress-energy 
tensor has the dimension of an energy density. 
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which allows the natural definition of the spin four- vector S A in such a way that 

ST = pL= e pvpa -^e-S* (2.7) 

y/-g A m A c 

where e pvpa is the four-dimensional antisymmetric Levi-Civita symbol such that e 0123 = 1. 
For the spin vector S A itself, we choose a four-vector which is purely spatial in the particle's 
instantaneous rest frame, where u A = (1, 0), hence the components of S A are (0, S A ) in that 
frame. Therefore, in any frame, 3 

Sy A = 0. (2.8) 

As a consequence of the supplementary spin condition (12.61) . we easily verify that 
di^S^ S A U ) / dr A = hence the spin scalar is conserved along the trajectories: S A ' S A = const. 
Furthermore, we can check, using (I2.6P and also the Papapetrou law of motion (I2.5p . that the 
mass defined by m 2 A c 2 = —p A p A is indeed constant along the trajectories: m A = const. Fi- 
nally, the relation linking the four-momentum p A and the four- velocity u A is readily deduced 
from the contraction of (12. 4p with the four-momentum, which results in 

fi(pu) A + m\t?vr A = 7^ 2 STS X /u A Ri Xp , (2.9) 

where (pu) A = p A u A . Contracting further this relation with the four-velocity one deduces 
the expression of (pu) A and inserting it back into (12.91) yields the desired relation between 
p A and u A . 

Let us from now on focus our attention on spin-orbit interactions, which are linear in 
the spins, and therefore neglect all quadratic and higher corrections in the spins, say 0(S 2 ). 
Drastic simplifications of the formalism occur in the linear case. Since the right-hand-side 
(RHS) of Eq. (12.91) is quadratic in the spins, we find that the four-momentum is linked to 
the four-velocity by the simple proportionality relation 

p\ = m A cu\ + 0(S 2 ) . (2.10) 

Hence, Eq. ( 12. 6 p becomes 

S%'itf = 0(S 3 ). (2.11) 

On the other hand, the equation of evolution for the spin, also sometimes referred to as 
the precessional equation, follows immediately from the relationship (12. 4ft together with the 
law (12TT0|) as DS p A /dr A = 0(S 2 ), or equivalently 

DS A 

= 0{S 2 ) . (2.12) 

dr A 

This is simply the equation of parallel transport, which means that the spin vector S A 
remains constant in a freely falling frame, as could have been expected beforehand. Of 
course Eq. (12. 12[) preserves the norm of the spin vector, S A S A = const. 

When performing PN expansions it is necessary to use three-dimensional like expressions 
(instead of four-dimensional) for the stress-energy tensor. The field point is accordingly 



The alternative choice S^p A = is equivalent to S^u^ = modulo cubic terms in the spins Q(S ) (see 



below) which are neglected in the present paper. Such choices are also adopted in Refs. [28l. |23. 136U37 1 
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denoted by x = (ct, x), and similarly the source points are denoted y A = (ct, ya)- The 
particle trajectories are considered as functions of the coordinate time t = x°/c, say y^(^), 
and we introduce the ordinary (coordinate) velocity v A (t) = dy A /dt, also a function of 
coordinate time. Using Eq. (I2.1(jp we can write the monopolar part (12.21) of the stress- 
energy tensor as 

= T"" + 0(S 2 ), (2.13) 



M 



NS 



where ns^^ is just the standard piece appropriate to point masses without spins, which 
reads, in three-dimensional form, 



rp jiV 

NS 



E 



m A 



V>A 



5(x- y A ) 



(2.14) 



We have referred to this part of the stress-energy tensor as the "non-spin" contribution (NS) 
in spite of its implicit dependence on the spins through the metric tensor. Here S = 5^ is 
the three-dimensional Dirac function. Similarly, the spin part of the stress-energy tensor, 
Eq. (12.31) . can be re- written as 



s 



S 



U A 



5(x-y A ) 
V~9a . 



[2.15) 



where the spin tensor S A (t) is now considered to be a function of coordinate time, like 
for the ordinary velocity v^(t). The covariant derivative V p acts on x, which appears in 
the argument of the delta- function as shown in (12.151) . and on time t through the time- 
dependence of the positions y A (t), velocities v^(t) and spins S^it). It is easy to further 
obtain the more explicit expression 



-g 

s 



A 



sfvp(x-y A ) +SZr% A v A 5(x-y A 



p(M -pu)A , 



(2.16) 



where Y vA = Y v pa (y A ) denotes the Christoffel symbol evaluated at the source point A, and 
where one should notice that the square-root of the determinant -J—g in the left-hand-side 
(LHS) is to be evaluated at the field point (t, x), contrarily to the factor l/y/— g A in the 
RHS of Eq. (12.151) which is to be computed at the source point y A = (ct,y A ). The explicit 
form (12.161) of the spin stress-energy tensor is used in all our practical calculations. 

In terms of three-dimensional variables the spin tensor reads [after taking into account 
the spin condition (12. 8p . namely Sq = —Sfv A /c] 



°A 



1 



„.A oA 

V~9A 



Jjk 



V^g~A 



oA , A^A oA 



where e^ k is the ordinary Levi-Civita symbol such that e 123 = 1. Here, we have 



(2.17a) 
(2.17b) 



U , 



fi'oo + 9oi 



(2.18a) 
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with u°a 



^jo + 9jk 



9pa V A Va/ c2 



(2.18b) 
(2.18c) 



In principle we could adopt as the basic spin variable the covariant vector (or covector) Sf. 
However, we shall instead use systematically the contravariant components of the vector 
S l A , which are obtained by raising the index on S£ by means of the spatial metric 7^, which 
denotes the inverse of the covariant spatial metric evaluated at point A, 7^ = g^- {i.e. such 
that YAlkj = dj)- Hence we define (and systematically use in all our computations) 



qi — ^ikqA 

D A — lA^k 



(2.19) 



Beware of the fact that the latter definition of the contravariant spin variable S\ differs from 
the possible alternative choice q A S^. The spin vector S A as defined by (I2.19P agrees with 
the choice already made in Refs. [36|, |37 . 



III. POST-NEWTONIAN METRIC AND EQUATIONS OF MOTION 



The starting point is the general formulation, i.e. valid for any matter stress-energy 
tensor T^ u with spatially compact support, of the PN metric and equations of motion at 
2.5PN order, as worked out in Ref. 

S3- 

In harmonic (or De Donder) coordinates 4 the 2.5PN 
metric is expressed in terms of certain "elementary" potentials as 



#00 = 

9oi = 

1 

9ij "ij 



1 H V V z H 

c 2 c 4 c 6 



V 3 

X + ViVt + — 

D 



(3.1a) 
(3.1b) 
(3.1c) 



These potentials, V, V$, • • • , are defined by some retarded integrals of appropriate PN 
iterated sources. To define them it is convenient to introduce the matter source densities 



a 



OA 



yuu _|_ ip 
( 

rpOi 

? 

C 

rpij 



(3.2a) 

(3.2b) 
(3.2c) 



(with T kk = 5ijT^). Then, with C^ 1 denoting the usual flat space-time retarded operator, 
we have for the Newtonian like potential V, 



V = □ r 1 {-4ttGct} =G J 



d 3 x' 



■ a fx', t — |x 



X' c 



(3.3a) 



Thus d u (v / ^5.9 A " y ) = 0) where g^ v is the inverse of the usu£l1 covaxicLirt metric Q^v-) 

and g = det(g p<7 ). 
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The higher-order PN potentials read 



(3.3b) 
(3.3c) 

(3.3d) 
(3.3e) 

All these potentials are subject, up to the required PN order, to the differential identities 



Vi 






% 


= D R 1 {-4 7 rGK- 


- 5ijcr kk ) - diVdjV} , 


Ri 


= □r 1 {-4ttG(Fct 


i - V { a) - 2d k VdiV k - 


X 


= Dp 1 {-AnGVa, 


+ 2Vid t diV + vdjv 




+ \{d t vf - 


2o,\ + ir^vjr} 









~W U + 2V 2 



dtV + dj 



Ri + VV 



oil 



W i:j - -6M 



ij vv kk 



(3.4a) 
(3.4b) 



which are consequences of the harmonic coordinate conditions; see Ref. [4j 

In this paper we shall specialize the latter PN metric to systems of particles with spin. 
In this case, as we have reviewed in Sec. [TTJ the stress tensor is the sum of the non-spin piece 
given by (12TT4D and of the spin part Q2H5J), thus T» u = NS T^ + S T^. Henceforth we often 
do not indicate the neglected 0(S 2 ) terms. Hence, the source densities (13.21) will be of the 
form a^ u = ns0>k + s<V>> an< ^ a ^ the potentials will thus admit similar decompositions, say 



V 



Wi. 



: v + v, ■■■ 

NS S 
NS J S J 



(3.5a) 
(3.5b) 



The equations of motion of spinning particles are obtained from the covariant conservation 
of the total stress-energy tensor, 







V U T^ + V u T^ + 0{S 2 ) 

NS S 



(3.6) 



To get the acceleration of the A-th particle, we insert into the conservation law ( 13. 6ft the 
expressions (12.141) and (I2.15P of the stress tensor, integrate over a small volume surrounding 
the particle A (excluding the other particles B), and use the properties of the Dirac delta- 
function. More precisely, in order to handle the delta-function, we systematically apply the 
rules appropriate to Hadamard's partie finie regularization and given by Eq. (14.61) below. As 
a result we obtain the equations of motion of the particle A and find useful to write them 
in the form 



dt 



(3.7) 



where both the "linear momentum density" P A and "force density" F A (per unit mass) 
involve a non-spin piece (NS) and the spin part (S), 



P: 



P A + P A 

N> s" 



(3.8a) 
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and F A 



F A + F A . 

NS P S P 



(3.8b) 



The non-spin parts correspond to the geodesic equations and read 



NS' 



NS f 



v v A g A 



nil 



1 V A V A {%9v\)a 

2 yf-gfrW* 



(3.9) 
(3.10) 



Their complete expressions in terms of the elementary potentials (13.31) were given in Ref. 47 
We shall need them for a spatial index (ji = i) and for completeness we report here the result 
(see Eqs. (8.3) in 0) 



NS 



-AVi + 3VV + -vV 



9, 



-8Ri + - VV + Wtjit - AVVi 
2 



+ -Vv 2 v* - 2v 2 Vi - AvVVj + 



+ o 



(3.11a) 



NS 



( + 



+ 



-VdiV + -v 2 diV - Av j diVj 
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AdiX + SVjdiVj - 8vidiRi + -v 2 VdiV 



7 



+ 2v j v k diW jk - 2v 2 v j d i V j + -v A dtV + -V 2 diV 



Av%diV - Av j VdiV 3 ] A + (Jzj . 



(3.11b) 



These expressions are still valid in the present situation, but we have to remember that the 
elementary potentials therein do involve contributions from the spins, e.g. V = ns^ + s^- 
Therefore it is crucial to compute the spin parts of the potentials and to insert them into 
the non-spin (geodesic-like) contributions to the equations of motion, Eqs. (13.111) . 

Now the purely spin parts, sP A and s-^f, will produce a deviation from the geodesic 
motion which is induced by the effect of spins. We have found that they admit the following 
expressions, 



- - n A V uA ^ p0 — -n A V 0A <J pv V ' A 
c^) pv 1 pa ^ A c^g^v 1 per J A „ ' 



m A c F A = - {d w g va ) A S P A U v\ 



-\V&*) A V ^ S A V°A 



(3.12a) 



(3.12b) 
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To compute them is relatively straightforward because all the metric coefficients and 
Christoffel symbols therein take their standard non-spin expressions (since we are looking 
for an effect linear in the spins), and these have already been computed in Ref. [47j ]. 

As a check of our calculations we have also used an alternative formulation of the equa- 
tions of motion, which is directly obtained from the Papapetrou equations of motion (12. 5p 
and reads, at linearized order in the spins, 

m A c^ = ~S x /u A R^ Xp + 0{S") . (3.13) 

We lower the free index /i so as to use the convenient relation Du^/cIta = du^/drA — 
^u A u\{d^g v \) a- The resulting equation takes the same form as Eq. (13.71) . 

dV A 

= (3.14) 

but with some distincts linear momentum and force densities V„ and It is clear that 
the non-spin parts, corresponding to geodesic motion, can be taken to be exactly the same 
as in our previous formulation, namely Eqs. ( 13. lip . However the spin parts are different; 
they are given in terms of the Riemann tensor -R^ Ao . T = R^xo^Ua) as follows, 

m A cV A = 0, (3.15a) 
s ^ 

J^A £ upar 

mAC & = ^ffro " A aL VA s " • (xl5b) 

where S A is the covariant spin covector appearing in (12.191) . The difference with Eqs. (5.1-3) 
in Tagoshi et al. j3?J is due to the fact that these authors work on the contravariant version 
of the Papapetrou equation. The advantage of the formulation (13.151) over the previous 
one (13.121) is of course that it is manifestly covariant. This advantage is however relatively 
minor in practical PN calculations, since the manifest covariance of the equations is anyway 
broken from the start. It remains that the two formulations are very useful, and their joint 
use provides a very good check of the calculations. 

The quantity (13.15bj) can be computed from the 2.5PN metric, by inserting it into the 



curvature tensor R^ XaT) but we may also express them directly by means of the elementary 
potentials (I3.3p . Let us give here the complete result at the required PN order, 

m A c F? = ^ {sijkfdjdtV + v l d jl V^jS k + 2£ jkl d u (Vv j - ^)^}^ 

'dfrV + v%v) (s k V + ^v 2 S k - {Sv)v k ^j + (d 2 t V + v l d i d t v)v i S h 

+ 2d l v[d l V k S j + S k d j V l + v j S k d l V^) + djVfdtV - v l d l V^jS k 
2{2d 3 Vd l V l - 2V j d il V + 20^8^ + Vd u Vj 




+ £jkl 



3 

2v j d l Vd l V + v j Vd u V + v j v m d lm V t - v j v m d u V m 



+ v^AVi - v'dAVj + v m duW jm - v m d lm W i3 - mWij - 2d il R j )S k 
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+ v 2 ( - duVj + ^d a v)s k + 2(Sv)v% l V j 



(3.16) 



IV. COMPUTATION OF THE SPIN PARTS OF ELEMENTARY POTENTIALS 

We shall compute all the spin parts of the elementary potentials listed in Eqs. (13.31) . which 
are needed for insertion into the "non-spin" parts of the momentum and force densities as 
defined by Eq. (13.111) . Here we do not compute the non-spin parts of the potentials since 



they are known from Ref. [47 



Let us start by deriving a few lowest-order results. First, it is immediate to see that the 
non-spin parts of the matter source densities a, Oi and o^, Eqs. ( 13. 2ft . start at Newtonian 
order, and that their spin parts start at 0.5PN order ~ 1/c in the cases of the vectorial 
erf and tensorial densities afj, and only at 1.5PN order ~ 1/c 3 in the case of the scalar 
density <r s . Here we are using our counting for the PN order of spins [see Eq. (Il.ip ]. which is 
physically appropriate to maximally rotating compact objects. With lowest-order precision 
the expressions of the source densities for two spinning particles read 

a^-^e ijk v i 1 S{d k S 1 + 1^2 + o(\) , (4.1a) 

o c \C J 

l - e ijk Si d k 6 1 + l++2 + , (4.1b) 



&i = 

s 2c 



S J 



-ie« (i v? SfoS! + 1 <- 2 + O (Jz) . (4.1c) 



The symbol 1 2 means adding the same terms but corresponding to the other particle. 
The Dirac delta-function is denoted by 5± = 5(x — yi), and dkS\ means the spatial gradient 
of S± with respect to the field point x. The lowest-order potentials are then straightforward 
to obtain from the fact that A(l/ri) = — Att 8\ (where r% = |x — yi|), and we get 

V = ~ s ijk v\ Si (I) + l 2 + (±) , (4.2a) 

Vi = -| Si d k (£)+l~2 + (D (1) , (4.2b) 

Wii = ~ e kl{i vf S* d l + 5ij e klm S{ d m (J^j + 1 <- 2 + O (±) , (4.2c) 

W kk = — e klm v\ S[ d m (- )+ 1 ^ 2 + O f ^ ) . (4.2d) 



s c V r i / V c 

At the dominant level only contribute to the potentials some compact-support terms (pro- 
portional to the source densities ozj) — notably the non-compact support term ~ dVdV in 

the spin part of the potential Wij, Eq. (13. 3d) , turns out to be negligible. 

To find all the spin terms in the equations of motion up to 2.5PN order, we see from 
Eq. (13.1 ip that we need V to 2.5PN order and Vi at 1.5PN order [i.e. 1PN beyond what 
is given by (14.2b|) ]. together with Wij, Ri and X at order 0.5PN. As we see, the potential 



Wj is already given by Eq. (I4.2cl) with the right precision. Our first problem is to obtain 
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the compact-support "Newtonian" potential V to the 2.5PN order. Definition (I3.3al) shows 
that the mass density a, source of V, admits at an arbitrary high PN order the structure 



a= /xi + fr )5 1 + -=dt( u l 5 l ) + -=d i [ i/^i +1^2. (4.3) 

The factors Jia-, s^a, s^a and ^u l A are functions of the spins and the velocities v A , and 
functionals of the metric components or, equivalently at 2.5PN, of the elementary poten- 
tials (13.31) . Note that though a^(x, t), by contrast to S A (t), depends on the field point, 
this is not the case of the moment-like quantities entering the square brackets of Eq. (12.151) . 
Each of them, being multiplied by the Dirac distributions 5a, is indeed evaluated at point 
x = yA, after the Hadamard procedure described below. Thus, it depends on time only 
(via the point-mass positions and velocities v^). The index S indicates an additional 
linear dependence in the spin components, but of course, the full spin dependence is more 
complicated due to the implicit occurrence of Sjf in the potentials themselves. Notably, the 
effective mass fii whose expression in terms of V, V*, Wu and v\ can be found in Ref. [47 
contains a net contribution due to the spin at the 2.5PN order and given by 

(ft) s ="^(-k + ?[- 4 r i " , '- 2 t- i ]) 1 +o (?) • (4 - 4) 

where the value at the particle's location is meant in the sense of Eq. (I4.6al) . The expressions 
of the other moments will not be provided here. It is in fact sufficient for our purpose to 
observe that, as shown by Eq. (14.41) . we have (/ii)s + s/^i — ^(l/c 5 ), an d that g^i is at least 
of order 0(1/ c 1 ) whereas $v\ is of order 0(l/c 3 ). 

As the spin contribution in a, say $0", is already of order 1.5PN ~ 1/c 3 , see Eq. (I4.1al) . 
we need to expand the retardations in V only at relative 1PN order, hence 

+ ^/^'|x-xl(^(x', t )) s + (^). (4.5) 

We then substitute the value of a following from Eq. (14.31) . The integrals are evaluated with 
the help of the formulas 

J d 3 x'F(x')5(x'-y 1 ) = (F) 1 , (4.6a) 

J d 3 x' F(x') %6(x! - yi ) = —(diF)i , (4.6b) 

where the values at point yi are denoted by parenthesis like for (F)\. These formulas extend 
the usual formulas of distribution theory, which are valid for a smooth function F with com- 
pact support, to singular functions with a finite number of singular points and deprived of 
essential singularities (see Ref. [48] for full explanations about this generalization). The for- 
mulas (14.61) are part of Hadamard's self-field regularization which is systematically employed 
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in the present approach and the one of Hdi. 5 In the end we are led to 



G 
G 



-Gui(t) U( 



— m 1 (a\) s d i r 1 + l 



y/-9{*>t) |x- x 



-A 



2 + - 



(4.7) 



The final result for V is obtained by replacing the moments and the determinant of the metric 
at the 2.5PN level by their explicit values derived from the lowest order approximation of 
the potentials. The computation of sVi is similar to that of sV, though slightly simpler since 
the counterpart of ft for cr, does not depend implicitly on the spin at the 1.5PN order. 

Next we explain how to compute the non-compact (NC) support terms, and we take the 
example of the particular NC term in the potential Ri given by 



R 

s 



(NC) 



A" 



-2d k Vd % V k 

s 



+ o 



(4.8) 



In the source of this term we have to insert the Newtonian approximation of the potential 



V, which is simply V 



+ 



Gm,2 



+ 0(c ), together with the leading-order spin term gVfc 



given previously in Eq. (|4.2bl) . The source being known we are then able to integrate (using 
the same techniques as in Ref. [47|) and we get 



J R (NC) 
s * 



G 2 m l 
8c 



G 2 m 2 



S 1 du d m g + 1 

1 2 



2 + 



in which 



satisfies 



g = In (n + r 2 + r l2 ) 



AS 



r x r 2 



(4.9) 



(4.10a) 



(4.10b) 



The crucial fact which enables the latter integration in closed analytic form is the existence 
of the function g (first introduced by Fock |5l|). This function and its generalizations are 
extremely useful in the computation of the spinless equations of motion at 2PN and 3PN 
orders 



Finally all the necessary spin parts of the potentials are computed by PN iteration, ready 
for insertion into the non-spin contribution of the equations of motion as given by Eqs. 03.1111 
For all the potentials we are in agreement with the results reported by Tagoshi et al. 
their Appendix. 6 



5 Hadamard's regularization is known to yield some ambiguous coefficients in the equations of motion and 
the radiation field of non-spinning point particles at 3PN order. When using dimensional regularization 
these ambiguities are seen to be associated with the appearance of poles oc 1/e (or "cancelled" poles) in 



the dimension of space d = 3 + e 12 1. The PN order considered in the present paper is merely 1PN, since 



we are computing the 1PN correction to the leading spin-orbit effect. At this order there are no poles; 
therefore dimensional and Hadamard's rcgularizations are equivalent. 
6 We have however noticed the following misprints in Ref. [37(: in Eq. (Alh) for s-Rj, the third term in 
the first parenthesis of the first line should be +1TI2/ (r^s 2 ); in Eq. (Ali) for gX, the first term in the 
parenthesis following (1112^2) in the third line must be read — mz/ {r\2S 2 ). 
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V. THE 2.5PN EQUATIONS OF MOTION WITH SPIN-ORBIT EFFECTS 



A. Equations in a general frame 



In addition to the spin parts of the potentials computed in Sec. [TV] and inserted into 
Eqs. (13. lip , we add the required spin corrections to the geodesic motion as given by either 
the formulation of Eqs. fl 3 . X 2 1) or that of (I3.15p -( l3.16l) . The latter corrections are computed 
by inserting into them the non-spin parts of the potentials taken from [47| . We find that the 
two formulations [respectively given by (13.121) and (I3.15I) - (I3.16I) ] are equivalent and agree 
on the result. Finally the 2.5PN equations of motion with spin-orbit effects are obtained in 
the form 



A N + -^A 



1 1 

1PN H o Al.5PN H 7 

c 6 s c 4 



A 2 PN + A2PN 

ss 

_ Gm 2 



1 

+ -r 



A2.5PN + A2.5PN 

S 



(5.1) 

Here the Newtonian acceleration is A l N = — ^r 2 -^, an d we denote by A l N , A\ PN , A 2PN and 

A % 2 5PN the standard non-spin contributions (in harmonic coordinates) which are well-known, 
see Eqs. (8.4) 



m 



47 



and earlier works reviewed in [52|. In particular A 2i5PN represents 
the standard radiation reaction damping term. (For simplicity we henceforth suppress the 
subscript NS on non-spin- type contributions.) 

The leading-order spin effect is the 1.5PN spin-orbit term. For this term we recover the 



standard expression, known from Refs. 23, 24( and given in in the center-of-mass 

frame, and in 37J in a general frame. In the following we shall sometimes use some formulas 
relating the "mixed products" of three vectors in three dimensions, 



(Ux,U 2 ,U 3 )V 



{UU X ) U 2 x U 3 + (UU 2 ) U 3 x Ui + (UU 3 ) Ui x U 2 
(U, U 2 , U 3 ) Ui + (U u U, U 3 ) U 2 + (U u U 2 , U) u 3 , 



(5.2a) 
(5.2b) 



valid for any vectors U, Ui, U2, U3 (in 3 dimensions). Here the vectorial product of ordinary 
Euclidean vectors is indicated with the x symbol, for instance (Ui x U 2 ) J = e % ^ k U{lJ\] 
parenthesis denote the usual Euclidean scalar product, (UU\) = U l U{ = U • Ui; and the 
mixed product, or determinant between three vectors, is denoted (U%, U 2 , U 3 ) = Ui • (U 2 x 
U 3 ) = e ijk U[UiUl. This yields 



A1.5PN 

s 



12 



'{3 1 ,ni2,v 12 ) , JS 2 ,n 12 ,v 12 ) 
h 



mi 



+ 6(ni 2 i;i2) 



n i2 x S 2 

m 2 



m 2 

V12 x Si 



mi 



ni2 + 3(n 12 f 12) 



L Vl2 X S 2 

m 2 



n 12 x Si 

mi 



(5.3a) 



We use, whenever convenient, the notation v 12 = vi — v 2 for the relative velocity. 

The next-order spin correction is the spin-spin (SS) at 2PN order. We do not give this 
term since we are concerned here with spin-orbit effects which are linear in the spins. The 
SS term is quadratic in the spins, 0(S 2 ), and can be found in Refs. |23|, [24| and e.g. in 
Eq. (5.9) of [37(] . Now the 1PN correction to the spin-orbit effect, which is the aim of this 
paper and the work [33], reads 



A2.5PN 

s 



Grrin 



12 



ni2 



6(n 1 2,v 1 ,v 2 ) 



(ui-Si) (v 2 S 2 ) 



mi 



+ 



m 2 



14 



+ Vi 



mi 

{S 2 ,n 12 , V12) 
m 2 

_ 3 (Si,n 12 ,v 1 
mi 



15(n 12 v 2 ) 2 + Q(vi 2 v 2 ) + 
I5(n 12 v 2 ) 2 + 6(vi 2 v 2 ) + 



26^1 + 18^ 



ri2 
49 Gm 1 



20 



ri2 
Gm 2 



(ni2«i) + ("12^2) + 6(7112^1) 



2 r 12 r 12 

{Si,n 12 ,v 2 ) 



-6(n 12 M (g2,Wl2,Vl) + (g2,Wl2,t,2) 



+ v 2 



6(ni 2 ui) 



m 2 

(S'ljTiiajVu) 



m 2 



I2(m 2 tfi) 



mi 
6(ni 2 w 2 ) 



-4 



, (£^^2) 
mi 



m 2 



(52,ni2,Ui 2 ; 



mi 



+ 



HL12 X Vi 

1112 x v 2 
1112 x Si 



3(n 12 v 12 ) 
6(ni 2 ui 2 ) 



Mi) 
mi 

(V2S 2 ) 

m 2 



6(ni 2 fi 

m 2 

Gmi (7112^2) 



+ 4 



4 



7*12 m 2 
Gmi (n 12 S 2 ) 
r 12 m 2 



Vi x v 2 



mi 



+ 4- 



m 2 



mi 
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y (7112^12) (7112^2) 2 + 3(7112^2) (v 12V2) 



14— (ni2«ia) 



7-12 



9— («12%) 
7-12 



+ 



+ 



+ 



ni2 x S 2 

m 2 

Gmi 
+ 

7-12 

V12 x Si 

mi 
V12 x S 2 

m 2 



15(7ii 2 wi 2 )(ni 2 f 2 ) 2 - 6(7ii 2 wi)(fi 2 w 2 ) + 12(n l2 v 2 )(v 12 v 2 ) 

Gm 2 



- y(ni2fi) + y(77i 2 7; 2 ) 



16- 



7-12 



1 (7712^12) 



15 G 

- 3(ni 2 i;i) (7712^2) + —(n 12 v 2 ) 2 H (14mi + 9m 2 ) + 3(t>i 2 7; 2 ) 

2 n 2 

„. , 9 w . 23 Gmi ,„Gwi2 
6(77i 2 t;2) + 4(t;i 2 7; 2 ) + — + 12 - 



2 r 12 



7-12 



(5.3b) 



Surprisingly, we find that our expression has substantial differences with the result given 
in Eq. (5.10) of 37J] - However, since we recovered in the last section exactly the same 
potentials as given in the Appendix of [13], and since as we shall see below we find perfect 
agreement with the equations of motion computed in [37| in the case of the center-of-mass 
frame, we believe that the latter differences can only be due to some trivial misprints (and 
most probably to some mixup of Mathematica files) in the last stage of the work 37] . 7 

In Appendix |A] we shall prove that the equations of motion stay invariant under global 
Poincare transformations. Such a verification is quite important to test the correctness of 
the equations (it played an important role during the computation of the 3PN non-spin 
terms in [1^, 5(|). Furthermore, we show in Appendix [B] that the test mass limit of the 



7 For completeness we indicate here all the misprints in Eq. (5.10) of [37]: in order to recover the 
correct acceleration, the last term before the closing curly brackets on the fifth line of Eq. (5.10), 
— \2e^ kl vP v\v l 2 (v\S\) /mi, must be replaced by — &e^ kl n^ v\v l 2 [{viSi) / rrii + {v 2 S 2 ) / rn 2 ]] the term before the 
last one in the seventh line has to be read — ftnPv^ '(nvi 2 )(v 2 S 2 )/m 2 instead of — %n?v\ (nvi 2 )(viSi) /mi; 
and the very last term +7v{v 2 (viSi)/nii must be modified as +v{v 2 [3(viSi) /mi + 4(i>2/S l 2)A n 2]. 
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equations of motion is identical with the geodesic equations around a Kerr black hole (for 
simplicity we restrict ourself to circular orbits). Both verifications have already been made 
in Ref. [37j but we present some alternative ways to do the checks. 



B. Equations in the center- of- mass frame 



Let us now present the result in the center-of-mass (CM) frame, defined by the nullity 
of the center-of-mass vector, equal to the conserved integral associated with the boost in- 
variance of the equations of motion, which will be checked in Appendix |A] We shall derive 
the center-of-mass integral at the 2.5PN order in the next section, however for the present 
computation we need it only at the 1.5PN order. When working in the CM frame we find 
it convenient to introduce the same spin variables as chosen by Kidder 29| (except that we 
denote by XI what he calls A), namely 



S 



(5.4a) 
(5.4b) 



Si + S 2 , 

/ S 2 Si \ 
ml ) . 

\m2 171%/ 

Mass parameters are denoted by m = ni\ + 7712, 5m = mi — m2 and v = m\ mj/m 2 (such 
that < v < 1/4). At the leading order in the spins we have the following relation between 
the positions yi and y2 in the CM frame and the relative position x = yi — y2 and velocity 
v = dx/dt = V! — v 2 = V12 (see e.g. Ref. [371 ]) 



yi 



Y2 



m 2 v 5m 

m 2c 2 m 

mi v 5m 

m 2c 2 m 



v 2 - 



Gm 



v 2 - 



r 

Gm 
r 



x + 



V 



m c 



3 vxS, 



x + 



v x S 



mc° 



(5.5a) 
(5.5b) 



In addition to the spin-orbit effect at order 1.5PN ~ 1/c 3 (last term in these relations), we 
have included the well-known 1PN ~ 1/c 2 non-spin term. This term is obviously needed 
here because during the reduction of the equations of motion to the CM frame at order 
2.5PN in the spins, we shall need to take into account the 1PN non-spin term coupled to 
the lowest-order 1.5PN spin term; such coupling evidently produces some 2.5PN spin terms. 
In the CM frame the equation of the relative motion reads 



d\ 1 1 1 

— — Bn + ^-Bipn + — B1.5PN + — 



c 3 s 



B2PN + B2PN 

ss 



1 



B2.5PN + B2.5PN 

s 



(5.6) 

where we recognize all the various terms similarly to Eq. (I5.ip . We find that the spin-orbit 
term and the 1PN correction to the spin-orbit are given, in terms of the spin variables (15.41) . 
by 



G t 

B1.5PN — —z\ n 



5 TTI 

12(S,n,v) + 6 — (E,n,v) 
m 



5m 

- 7v x S - 3 — v x E 

m 



G , 

B2.5PN — —z \ n 



5tyi 

+ 9(nv)n x S + 3 — (nv)n x S 
m 



Gm 



(5.7a) 



(S, n,v)[ - 30u(nv) 2 + 24vv 2 —(38 + 25u) 



16 



H (S,n,f) - lbv(nv) 2 + 12z/w 2 (18 H v 

m \ r 2 



+ (m;)v 
+ n x v 



(S,n,v)(-9 + 9u) + —(E,n,v)(-3 + 6v) 



Gin 



Sm ( . Gm 



(nv)(vS)(-3 + 3v) - 8 u(nS) - 

r m \ r 



i/(nS) + 3{nv)(vY) 



+ (nv )n x S 



_ 1^(^)2 + 2 W - ^(25 + 15z/) 
2 r 



5m . . 

■ (nv )n x S 

m 



- lSz/fnt;) 2 + 12z/f 2 



Gm . 17 , 

9 + -Z/ 

r 2 



+ v x S 



2 r 



<5m 
■ v x S 

m 



9z/(ra;) 2 - 7vv 2 + —(9 + -v) 
r 2 



(5.7b) 



We find perfect agreement with Eqs. (5.18) and (5.20) of Tagoshi et al. [3 



C. Reduction to quasi circular orbits 



Finally we present the case where the orbit is nearly circular, i.e. whose radius is constant 
apart from small perturbations induced by the spins (as usual we neglect the gravitational 
radiation damping at 2.5PN order). Following Ref. [29j] we introduce an orthonormal triad 
{n, A, £} defined by n = x/r as before, £ = Ln/|Ln| where Ln = px x v denotes the 
Newtonian angular momentum, and A = £ x n. The orbital frequency u is defined for 
general, not necessarily circular orbits, v = fn + ruo\ where f = (nv). The components of 
the acceleration a = dv /dt along the basis {n, A, £} are then given by 



n • a 
A a 



ruo 



dt) 



(5.8a) 
(5.8b) 

(5.8c) 



We project out the spins on this orthonormal basis, defining S = S n n + S\X + Sg£ and 
similarly for S. Next we impose the restriction to circular orbits which means r = = r 



and v 2 



2, ,2 



r uj 



'neglecting radiation reaction damping terms). In this way we find that the 



equations of motion (15. 6p with (15.71) are of the type 



dv 9 

— = — urr n + 

dt 



(5.9) 



Note that the spin variables adopted in [37( are defined by Xi 
and differ from our own. We have 



_ l 

— 2 



(1 



„ . Sm 

^v)Xs H Xa 

m 



and £ 



Si 



dm 



and Xa 



( Si _ s 2 A 

\™1 ™2/ 



17 



There is no component of the acceleration along A. Comparing with Eqs. (15.81) in the case 
of circular orbits, we see that u is indeed the orbital frequency, while ae = —ruo{\ ■ di/dt) is 
proportional to the variation of £ in the direction of the velocity v = ru\. We find that uo 2 
is of the form 



Gin 



1PN 



1 1 

-Z Cl.5PN + ~T 



C 3 s 



C2PN + C 2PN 
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C2.5PN ( 




+ lf 


♦•(i) 


ss 
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(5.10) 



where Cipn and C2PN denote the standard non-spin contributions, 9 and where 



G.5PN 
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C2.5PN 
s 



Gm 
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39 

y 



8m _ 
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m 
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v)St + 



21 11 

Z7 
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5m 



m 



(5.11a) 
(5.11b) 



with, e.g. Sg = (S£) = S • I. On the other hand, we get 



1 1 1 

«1.5PN + -7 « 2PN + -7 «2.5PN + C 



C 3 S 



C 4 SS" 



with spin-orbit coefficients 



c^i.spn 
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(5.12) 



5/2 



mr z 



63 zA 27 5m 

"T + 2 J 6n ~Tm n 



(5.13a) 
(5.13b) 



We see that the resulting motion cannot be exactly circular for general orientations of 
the spins. Let us show however that the time-averaged acceleration coincides with the accel- 
eration of a particle that rotates uniformly about the origin. In a first step, we must make 
explicit the time dependence of the dynamical variables x, S and S. As the motion is uni- 
formly circular in the absence of spin, the position x decomposed along a fixed orthonormal 
basis {e 1 ,e 2 ,£} reads 

x(i) = eir cos^ns^) + e 2 r sin(a;Ns^) , (5-14) 

with cg>ns being the orbital frequency when the spins are turned off. 

The spin variables are computed by means of the precession equations, which decouple in 
the case of a pure spin-orbit interaction. The spin 1, for instance, obeys an equation whose 



right-hand-side is polynomial in Gm/r 



(nSi) and (vSi). For dimensional reasons, it 



must then have the form (for circular orbits, up to say the 2PN order) 



dS\ 
~dt 



fc=l,2 



/Gm\ k 
\ rc' 



[ rr 2 ) 



(k,n) 



vSi)n + a^ v \nS\)v 



+ 



(5.15) 



and similarly for d^/dt. The functions of mi/m, mj/m denoted by Og and a^ v ' > may 
be obtained from the results of the next section (see also paper II). They allow us to define 



They are given by Cipn 



Gm 



(—3 + v) and C2PN = (^r 1 ) (6 + ^rV + v 2 ) in harmonic coordinates. 
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dimensionless coefficients like a s " = Ylk=i 2(Gm/r) k a! y ^ n \ The key point is that the latter 
coefficients are constant, which suggests to solve the above differential equations in the 
moving basis {n, A, £}. Indeed, the time derivative of a spin component in this basis, say 
= (nSi), is given by a relation of the type 

= (nSt) + (nSd (5.16) 

at 

with ri = x/r = c^ns-^- This results, after eliminating Si by means of Eq. (15.151) . in a linear 
differential equation with constant coefficients for S„. Proceeding in the same way for the 
other components of the first spin, we arrive at the following system: 

-^A = M 5l .X 5l , (5.17) 

where is a 3 x 3 constant matrix and = (S^,S\,Sj). The relations £.dS\/dt = 
E.dSz/dt = (since £ is constant because we neglect the SS terms) imply that (0,0, 1) is 
an eigenvector associated with the eigenvalue Ao = 0. There remain two eigenvalues, say 
A^ and A^; but since the trace of M Sl vanishes because (nv) = 0, we have \± = — A^. 
In the end, we notice that the spins are almost constant at Newtonian order in the basis 
{ei, e2, £}, which means that they precess about £ with angular velocity — c^ns i n the moving 
frame. Therefore, Xf is purely imaginary and reduces to ±2ix>ns at Newtonian order. At 
higher order we shall have Xf = ±«(u;ns — ^i) where fii = 0(1/ c 2 ) represents the precession 
frequency. The components and S\ solving Eq. (15.171) are then linear combinations of 
cos[(— cj NS + f2i)t] and sin[(— u;ns + ^i)*]- 10 As for the component Sj, it is constant neglecting 
terms quadratic in the spins. 

We complete our proof by time averaging the term ae in the acceleration (I5.9p . We 
first observe that the conservative part of the dynamics involves three different angular 
frequencies (c^ns, &i and f^), so that it cannot be periodic in general. Therefore, it is not 
appropriate to average the particle motion on the orbital period. Instead, the time-average 
will be achieved on infinite time. Defining 

i r t+T 

(Si) = Jim - / dt' Sl{t'), (5.18) 



we find (S^) = 0. We next notice that the orbital frequency uj is actually constant (neglecting 
SS terms), for it depends on the spin through Sj and Sj only, which are constant. The 
average of at is a linear combination of (S n ) = ((nS)) = and (E n ) = = 0; hence it 

does not contribute: (at) = 0. 



VI. THE 2PN SPIN-ORBIT EQUATIONS OF PRECESSION 

In this Section we give the equations of evolution of the spins, or precession equations, at 
relative 2PN order, i.e. one PN order beyond the dominant term. The precession equations 



This can also be deduced immediately from introducing a different spin variable SI with constant magni- 
tude (described in Sec. VII of paper II) and obeying dSl/dt = £1% x S^; noticing that the components of 
S^j; in the basis {n, A,^} are linear combinations of those of Si, with constant coefficients. 
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are quite simple to derive from the equation of parallel transport (I2.12p . which we recall is 
valid at the linear order in the spins [neglecting O (S 2 )], but at any PN order in that term 
which is linear in the spins. The PN corrections are easil y co mputed from the non-spin part 
of the metric and Christoffel symbols computed in Ref. [47]] . The precession equations in a 
general frame take the form 



dSi 
~dt 



1 



1 



~Z TlPN H n T1.5PN H 7 T2PN + O — 



c 2 s 



c 3 ss 



c 4 s 



(6.1) 



together with the equation with 1 

Grri2 



TlPN 

s 



12 



2. At the lowest order we find 
Si(ni 2 t>i 2 ) - 2n 12 (v 12 Si) + (vi - 2v 2 )(n 12 S 1 ) 



(6.2) 



The above equation is already known [29|, [37j. See e.g. Eq. (4.3) in 37] and the paragraph 
afterward commenting about the difference with formulations based on an alternative defi- 
nition for the spin, like that of Ref. j29[. The spin-spin (SS) term is also known but is out 
of the scope of the present paper (and the parallel transport equation we employ); it can be 
found elsewhere, see Eqs. (2)-(3) of [8|. Then we find that the next-order spin-orbit term is 



T2PN 
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12 
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G 
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ri2 
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v 2 



+ 2(n 12 v 12 ) ( (wiS 1 !) + (^^i) 



-{ni2Si){n 12 v 2 ) 2 + (vi 2 Si)(ni 2 v 2 ) 



{nx 2 S x ) [2(v 12 v 2 ) + 3(n 12 v 2 ) 2 
G 

5(ni 2 5i) — (mi - m 2 ) 
r\2 



(6.3) 



For completeness and for the benefit of users of these formulas in the data analysis of 
detectors, we present also the precession equations in the CM frame, using our specific spin 
variables defined by (15.41) . These are 



dS 
~dt 

~dt 



11 1 
Uipn H — o U1.5PN H — 7 U2PN + O 



c 2 s 



c 3 ss 



c 4 s 



11 1 

Vipn H — ? V1.5PN H — j V2PN + O 



c 2 s 



c J ss c 1 S 



1 



(6.4a) 
(6.4b) 



where all the spin-orbit coefficients are given by 
G m v 



UlPN 
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U2PN 
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-A(v3) - 2— (vE) 
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+ v 
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(6.5a) 
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and 
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(6.6b) 



To these expressions one may add the SS terms in the standard way (see Eqs. (2)-(3) of pi). 



VII. SPIN EFFECTS IN THE CONSERVED INTEGRALS OF THE MOTION 

Having obtained in Sec. [V] the equations of motion, the important task is now to deduce 
from them the complete set of conserved integrals of the motion associated with the global 



Poincare invariance of these equations (which has been checked in Ref. 37]] and Appendix lAl 
below). In principle, the conserved integrals of the motion, which generalize the usual 
notions of energy, angular and linear momenta, and center of mass position, should be best 
derived from a Lagrangian. In the present paper, however, we did not attempt to derive 
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a complete Lagrangian for the particles with spins (see 43j for a discussion on how to 
formulate Lagrangians with spins); rather, we have obtained the integrals of the motion 
by "guess-work" , starting from their most general admissible form, and then imposing the 
conservation laws when the equations of motion are satisfied. 11 Here we simply state the 
results. 

The PN expansion of the conserved integral of the energy, namely E such that dE/dt = 0, 
reads as 



E — Eqy H — j-E-ipn H — n Ei.spn H — r 



c 3 s 



-^2PN + -E2PN 
SS 



+ -T E2.5PN + 0[- 



c 5 s 



(7.1) 



where the non-spin pieces, Ej$, E 1PN and -E2PN, are known and can be found e.g. in Ref. [54 



For instance we have Ej$ = \vrt\ v\ + \rri2 v\ — G ™^™ 2 . For the lowest-order spin-orbit effect 



we find, in agreement with the standard result 

Gm,2 
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(Si, n 12 ,ui) + 1^2 



(7.2a) 
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where we employ a special notation for the totally anti-symmetric "mixed product" between 
three vectors, as given in (15.21) . For the spin-orbit contribution at 2.5PN order we find 
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ri2 

, Gm 2 
ri2 

(7.2b) 



Notice that several equivalent forms can be given to this result. For instance if wished 
one could introduce the mixed product (Si, f 1,1*2) in place of a (ni 2 , 1*1, i* 2 ) in the last term 
of ( 17.2bl) . making use of linear combinations such as (ni 2 i*i)(Si, 1*1, 1*2) = (7^12, 1*1, 1*2) (1*1 Si) + 
(Si, ni2, 1*2)1*1 — (Si, ni2, vi)(viv 2 ) [a consequence of Eq. (15. 2j) ]. As before we do not give the 
SS contribution at 2PN order (see [29] for instance). 

We give here the corresponding result for the conserved center-of-mass energy in the CM 
frame: 



E = mu c' 



r 1 1 1 

|e N + — eip N + — ei.5PN + — 



C2PN + g ^2PN 



+ ^e 2 . 5 PN + O 



C° S 



(7.3) 



where eN 



12 _ Cm 



see [55| for the other non-spin contributions). The SO coupling terms 



in the CM frame are found to be 



£l.5PN : 

s r 



— <(S,n,v) + (E,n,i*^ — 



m 



[7Aa) 



As usual we neglect the radiation reaction effect at 2.5PN order. Indeed we know that such an effect does 
not depend on the spins. The contribution of the spins in the radiation reaction force comes in at 1.5PN 



order beyond the dominant effect, which means at the 4PN level, and has been computed in Ref. [53 
Radiation reaction effects will be included into the present formalism when we obtain the contributions 
of the spins in the GW flux [38 1. 
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Let us next deal with the conserved total angular momentum J, i.e. dJ/dt = 0, sum of 
orbital and spin contributions, which we write as 



1 1 
L + -Sx + -S 2 

c c 



(7.5) 



where L is the orbital angular momentum, and where Si and S 2 are the contravariant spin 
vectors defined following the specific choice made in Eq. (I2.19P [recall also their peculiar 
dimension which follows from (11. ip ]. The angular momentum L admits the PN expansion 



_ 1 _ 1 1 
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(7.6) 



where all the non-spin pieces are given by Eq. (4.4) of [54|. For instance, Ln = m\ yi x vj + 
m 2 y 2 x V2. Now, in order to express in the best way the spin-orbit contributions in L, we 
find that they must be written in the following way, 



L1.5PN = Yi x Pi + Y2 x P 2 + K1.5PN 

s s s s 

L2.5PN = yi x q x + y 2 x q 2 + K2.5PN • 

s s s s 



(7.7a) 
(7.7b) 



in which we have introduced some convenient notions of the "individual linear momenta" 
of the particles, say sPi and sP2 at 1.5PN order, and sqi and sq2 at 2.5PN order. The 
extra terms in the RHS, sKi.spn and sK 2 .spn, incorporate all what remains, the point being 
that they depend on the positions of the particles only through their relative separation, 
i.e. t\i = |yi — y 2 | and ni 2 = (yi — y 2 )/ri 2 . The only dependence of the conserved angular 
momentum on the individual positions yi and y 2 is the one which is given explicitly by the 
first terms of Eqs. (17.71) . 12 The results we find for these momenta are 
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+ 1112 x V12 



3(ni 2 S , i)((ni 2 Vi) + {n 12 v 2 ) ) + (i^Si) 



(7.8b) 



12 However, let us stress that the definition of some individual momenta for the particles is merely introduced 
here as a convenient notation. In order to define in a meaningful way the notions of individual linear 
momenta of the particles with spins, we would need a Lagrangian, which as said before we did not compute, 
and the linear momenta would simply be the conjugate momenta of the ordinary positions. 
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together with the equations with 1 «-> 2. The last terms in the RHS of Eqs. (17.71) are 
explicitly given by 
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The 1.5PN term in the conserved angular momentum, Eq. (I7.7aj) . agrees with the result of 



Kidder 29 



13 



Let us add a comment on the meaning of the conservation of the total angular momentum 
J at 2.5PN order [Eq. (17. 5p with (I7.6p ], When differentiating J with respect to time, we 
generate several spin contributions at 2.5PN order: (i) The "main" one is coming from the 
differentiation of the Newtonian term Ln, and is due to the replacement of the acceleration by 
the equations of motion (15.11) with (]5.3b|) ; (ii) There is the one coming from the differentiation 
of the 1PN part Ljpn, since the replacement of the accelerations at order 1.5PN [Eq. (I5.3al) 
therein does also produce some terms at 2.5PN order; (iii) When differentiating the lowest- 
order spin-orbit term sLqpn, the derivative of the spins gives other 2.5PN terms via the 
precessional equations; (iv) When differentiating the spin vectors themselves, Si and S 2 , 
one must make use of the precessional equations with their full 2PN accuracy 14 which are 
given by Eqs. (I6.ip - (I6.3I) . Only when account is taken of all these replacements (i)-(iv) of 
accelerations and spin precession, does one find that J is conserved, dJ/dt = 0, up to 2.5PN 
order (neglecting the 2.5PN non-spin radiation reaction damping). 

The orbital angular momentum in the CM frame reads: 
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(7.10) 



where £® = mx x v; the non-spin contributions can be found in Refs. [5J, |55|, [56] . We have 

• 2 Gm s 



1 6m 9 _ iv 
s 2 m \ 2 
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Ref. 29( uses different definitions for the spin variables, which are related to ours by 



rv2j 2c 2 



14 



This is the only place where one needs the precessional equations with 2PN accuracy. 
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Finally let us give the conserved integrals of the linear momentum P and center of mass 
position G, which are related to each other by dG/dt = P. Recall that the existence of 
the center-of-mass integral G is a consequence of the boost-invariance of the equations of 
motion (cf. Appendix |A|) . Both P and G admit a PN expansion exactly like those of E 
and L. Quite naturally, we find that the spin-orbit contributions in P are simply given by 
the sum of the "individual" linear momenta for each particles that we found convenient to 
introduce in order to express the angular momentum in Eqs. (17.71) . Thus, 



Pl.SPN = Pi +P 2 

s s s 
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s s s 
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where the explicit expressions (17. 8p hold. For G, we obtain rather simple expressions, 
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The derivation of the complete set of integrals of the motion gives us further confidence 
in the physical soundness of the equations of motion derived in this paper. Those results, 
together with the analyses performed in Appendices \M and [Bj complete the resolution of 
the problem of linear spin-orbit effects in the binary's equations of motion at 2.5PN order. 
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APPENDIX A: LORENTZ INVARIANCE OF THE EQUATIONS OF MOTION 



Because of the global Poincare invariance of the Einstein equations (with bounded 
sources), and the manifest covariance of the De Donder harmonicity condition, it is not 
possible to physically distinguish between two harmonic-coordinate grids differing by a mere 
Lorentz transformation. As a result, the equations of motion must be of the same form in 
two such grids. In other words, up to an arbitrary PN order n, the link between the boosted 
acceleration a'^yc, vc, a c) and the boosted positions y' B (yc, v c), velocities, v' B (yc, v<?) and 
spins S^yc, vc, Sc), must be given by the original equations of motion [i.e. Eq. (15. ip at 
the 2.5PN level] with the original variables being replaced by their primed counterparts. 
Note that the Euclidean metric and the totally antisymmetric tensors remain unchanged 
under Lorentz transformations. Schematically, we may write a[ = A(y' B ,v' B ,S' B ,5ij,Eijk) 
for B — 1,2. The resulting relation between un-boosted quantities, 

ai(yc,vc,a c ) = A^(y c , v c ), v^(y c ,v c ),S' i j(yc, v c , S c ), fc^ey*) +C>(^ T ) , (Al) 

defines a function A' as ai = A'(yc, vc, Sc, + 0{l/c n+1 ). Equivalence with the 

equations of motion in the un-boosted frame: ai = A(ys, Vb, &b, £ijk) + Oil/ c n+1 ), 
means precisely that 

A = A' , (A2) 

up to negligible PN corrections. This property constitutes the so-called explicit Lorentz 
boost invariance of the equations of motion. It happens to be a very powerful check for the 
coefficients entering the functions A of Eq. (15.11) . and in particular its contribution A 2 .5pn 
[see Eq. (IPb|) ]. 

In order to verify the validity of Eq. ( 1A2[) . we need to determine the function A' explicitly, 
which requires to know how y#, vg, ai and Sb transform under a Lorentz boost. Let us 
start with considering an arbitrary space-time event P with coordinates x M in the current 
working frame (J 7 ) . Its coordinates in a boosted frame (J 7 ') of relative velocity V are related 
to the original ones by x' M = N 1 l/ (V)x u , where the Lorentz matrix A^(V) is given by 

A° (V) = 7, (A3a) 
A (V) = A° l (V) = - 7 ^, (A3b) 

J 7 + 1 c z 



with 7 being the Lorentz factor 1/y 1 — V 2 /c 2 . An event Q with coordinates y 1 ^ in (J 7 ') is 
simultaneous to P in the new frame if and only if y'° = x'°. There exist two such events 
located on the two world-lines of the binary companions. Their coordinates in (J 7 ') are 
denoted by y'± = (ct',y[) and y'^ = (ct',-y' 2 ) respectively. The mapping t' — > y^ defines a 
function y[(t'), and similarly for the second body. The events having coordinates (ct', y'x(t')) 
and (ct', y-iit')) in (J 7 ') do not generally appear as simultaneous in (J 7 ). They may be referred 
to in components as (cti, yi(ti)) and (ct 2 , y2(^)) in that frame, the functions yi(t) and y 2 (^) 
being the original trajectories. By construction, we have 

y'm = AP„(V)tf ft) . (A4) 
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Let us express in the end the RHS in terms of the coordinate time t. A derivation of the 
general formula linking y, (f) to yi(t) in the PN scheme can be found in 57J. This relation 
reads, see Eqs. (3.20) in 



yi(0 = yi(*)-TV(t 
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c 2 7 + 1 
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The velocity and acceleration follow from the partial derivation with respect to t' together 
with the formula d' t = 7 St + r yV t di. 
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The spin components in the new frame cannot be obtained directly from the linear Lorentz 
transformation law. This is because the definition of Si and S2 involves the inverse of the 
3-metric 7^ induced by g av on a slice t = const. Now, 7^ implicitly depends on the choice 
of the coordinate time and is generally singular because of the particle's self-gravitation. 
To avoid complications rising from this second issue, we shall first focus on the case of test 
particles on a fixed background. 

In the frame {J 7 '), the spin components of the first test body read 



(A7) 



with Ji 3 (t') = 7'* J (y' 1 , t'). Whereas the transformation law of 7/ is more difficult, that of 
(7y)i = (9ij)i results straightforwardly from the transformation of the space-time metric: 



(^■)i(0=A^(V)A^(V)(^) 1 (t lj , 



(A8) 



with A A ^(V) = A M A (— V) denoting the inverse transformation. Therefore, computing 7 1 ? ' ? (i / ) 
amounts to expressing the latter quantity as a function of (5+)i- This is achieved by means 
of the relation det(7^)i(7' 4 - ? )i = (Conr/)/, valid for any matrix (7^)1 between its deter- 
minant det(7^)i, its comatrix (Conr/)^ 7 and its inverse. For 3-dimensional matrices, the 
determinant may be written in an Euclidean covariant form as 



det( 7 yi = 

Similarly, we have for the comatrix 
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The inverse spatial metrics 7^ is then given by the ratio of the RHS of Eqs. (IA10I) and (IA9I) . 
where the primed metric relates to (<? M i,)i after Eq. (1A81) . 

We finally look at the determination of the covariant spin components S*/. As is a 
Lorentzian vector, they are at once seen to be equal to 



sf(0 = V(v)si(fi) 



(All) 
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and, by virtue of the supplementary condition (12. 8p . Sq = —S}v\/c. 

At this stage, we have expressed S'\ in terms of quantities evaluated at time t±, which 
has led us to a relation of the form Si = <Si(ti). It remains to rewrite «Si(ti) as a function 
of t. For the present purpose, we restrict ourselves to a perturbative approach, and resort 
to the convenient formula 



l) n 



,71 — 1 



n=l 



(A12) 



generalizing in a straightforward way Eq. (3.16) of Ref. [57] to any smooth function / (see 
also the Appendix A of [H3|). In the end, this yields the following identity for the spin 
"vector" S[ = (S*) defined in the frame (J 7 '): 
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rfgl 
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where 
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These expressions are valid at any order in the boost velocity V. After specializing the 
above equation truncated at the PN level to the metric (13.11) . we arrive at 



s; = s x + ^ ( - fasi) + ^(v^)) + o(4 



(A15) 



Note that all powers of V consistent with the 1PN approximation beyond the leading 
spin-orbit term have been included. In principle, Eq. (lA14aj) holds only for test parti- 
cles. Nonetheless, it turns out not to depend on any regularized field. It is thus legitimate 
to extend it to the conditions of the present problem. 

With the previous transformation laws in hand, we are in position to check the Lorentz 
invariance as explained before. After a lengthy calculation, we arrive at the expected identity 
A(yBj v b> Sb, 5ij, Eijk) — A (yb, ye, Sb, e^k) = 0. 



APPENDIX B: TEST-MASS LIMIT OF THE EQUATIONS OF MOTION 

In the limit where one of the objects, say the number 1, is nearly at rest while its 
companion has a very small mass for a finite ratio S2 /tti 2 , we must recover the dynamics 
of a spinning test particle in the background of a Kerr black hole of mass m\ and spin 
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Si = m\a\ (in this Appendix we pose G = c = 1). To allow direct comparison with the PN 
equations of motion for m-i — > at S2/7712 = const, we shall work with the Kerr metric in 
harmonic coordinates. The link between the Boyer-Lindquist grid (indicated by the label 
BL henceforth) and some spatial harmonic coordinates can be obtained from Eqs. (41) and 
(43) of Ref. Q: 



x 1 + ix 2 



x 



^bl — mi + ia\ ) sin 6bl exp i 
(t B l - mil cos^bl , 



+ 



a 1 



In 



r + — r_ 



tbl - r. 



(Bla) 
(Bib) 



with r± = mi ± ymf — a\ and i = —1. Since V^V^tjbl = 0, we may also choose t = £bl- 
The exact expression of the metric in the new grid is rather complicated, but we shall not 
need it beyond the linear order in the spin. Neglecting the quadratic terms O(Sf), the line 
element reduces to 



ds 2 



mi 2 Amiai . 2 r + m x 2 
at sin vdtdcp H ar 



- 2 



mfai r + m\ 
r 2 r — rrii 



r + mi 



r — mi 

2 / j/i2 



sin 6drd(f) + (r + nil) (d0 + sin 



(B2) 



which coincides with the one deriving from the metric (13.1 1) at the dominant order, hence 
the harmonic coordinates defined by Eqs. fIBlj) and t = £bl are the same as those of the PN 
formalism. 

At this level, we may derive the equations of motion of a test particle with spin per unit 
mass 5*2/777.2 orbiting in the gravitational field flB2j) . For simplicity, we assume the trajectory 
to be circular and lie in the equatorial plane 9 = 7r/2; the vector d z points to the direction 
of the spin black hole, so that Si = Sf = niiai; the spherical coordinate basis is denoted by 
(d r , de, 8$). The circularity conditions state in particular that r remains constant in time. 
The spatial components of the four-velocity are then 
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After taking these relations into account, the explicit form of the evolution equations ( 13. 14[) 
becomes 
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d 

dt 



u 



g<t>o + 



dt 







(B4d) 



The harmonic gravitational field only depends on r and 9, both of which do not change with 
time. It is itself independent of t. Thus, Eqs. flB4aj) and flB4dj) imply that u° and to = d<p/dt 
are constant, whereas flB4cl) yields S 2 = 0; flB4bl) shows that Sf = const and fixes the 
value of uj. We draw the time variation of the spin S2 from the precession equation (I2.12p 
specialized to the Kerr background (IB2I) : 
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Noticing that dS 2 /dr = 0, it is immediate to see from ( 1B5al) together with the condition 
SqU° = —S^u^ that S| = 0. The remaining equations are identically satisfied. As a result, 
the spin of the small object is aligned (or anti-aligned) with the spin of the black hole, 
meaning that 



Side 



;S 2 d 7 



(B6) 



(r + nil) 2 

up to possible quadratic contributions. In the test particle limit, the spin vectors are related 
to S and S as Si = S + C(m 2 ) and S 2 /m 2 = (S + S)/m + 0(m 2 ). Insertion of these values 
in Eq. ()B4b|) leads to the solution 



9 m 
uj = — 
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7 



3/2 



r 3 Ul + 7) 3 m 2 (l + 7) 9 / 2 



5S Z + 3S Z + 3 7 (S 2 + S z ) + 0{S 



(B7) 



with 7 = m/r = nii/r + 0{m 2 ). By expanding the latter equality at the 2.5PN order, we 
recover the generalized Kepler relation given by fl5.10|) - fl5.11l) for v — > 0. 
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